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Coulomb blockade effects in capacitively coupled quantum 
dots can be utilized for constructing an A'^-qubit system with 
antiferromagnetic Ising interactions. Starting from the tun- 
neling Hamiltonian, we theoretically show that the Hamil- 
tonian for a weakly coupled quantum-dot array is reduced 
to that for nuclear magnetic resonance (NMR) spectroscopy. 
Quantum operations are carried out by applying only electri- 
cal pulse sequences. Thus various error-correction methods 
developed in NMR spectroscopy and NMR quantum comput- 
ers are applicable without using magnetic fields. A possible 
measurement scheme in an A'^-qubit system is quantitatively 
discussed. 

03.67.Lx,73.23.-b 



Quantum computers have been widely investigated 
from many perspectives Q-O]. A quantum-dot array 
is a promising candidate for the basic element of a quan- 
tum computer from the viewpoint of technological feasi- 
bility |]1^,|4| . It is noteworthy that controlling of lateral 
and vertical order in self-organized quantum-dot super- 

Although a spin-based 



lattices has been realized |15 



quantum dot computer has been intensively discussed 
|g|,|| , it seems that the quantum dot system using charged 
states is more accessible because the latter will be able 
to be constructed of various materials other than III- 
V group materials such as GaAs. From this perspec- 
tive, several authors investigated the quantum computer 
based on the charged states ]10|-[l^. In Ref. ||l^, we 
discussed the quantum dot computer in the limit of a 
free-electron approximation. The free-electron approxi- 
mation will be valid only when the interdot tunneling is 
strong and Coulomb blockade is suppressed, or the size of 
the quantum dots is as small as that of artificial atoms. 
Generally speaking, however, given the current state of 
technology, it seems that it will not be possible to make 
quantum dots as small as atomic order in the near fu- 
ture. Moreover, controUed-NOT (CNOT) operation for 
a quantum computer based on charged states has been 
discussed only in the two neighboring qubits so far; a 
scheme for constructing more than two qubits on a de- 
vice remains unclear. 

In this paper we advance the analysis of the quantum 
computer based on charged states and show a general 
A^-qubit scheme. Starting from tunneling Hamiltonian 
in the Coulomb blockade regime, we demonstrate that 
the Hamiltonians for one- and two-dimensional arrays 
of weakly coupled quantum dots are reduced to those 
for standard nuclear magnetic resonance (NMR) spec- 
troscopy p7[| . This enables any quantum computation to 



be described by electric pulse sequences. We assume that 
interdot tunneling is not so strong that the capacitance 
between two quantum dots is defined and we can use the 
Hamiltonian which allows the electrostatic energy of the 
quantum-dot structure to be described in terms of the 
capacitances of individual dots and gate voltages. This 
type of quantum-dot system was discussed experimen- 
tally by Livermore et al [|6|. The interaction between 
qubits is derived from the Coulomb interaction between 
the excess electrons in the quantum dots. We also nu- 
merically illustrate the measurement process using the 
field-effect, which is considered to be a sensitive mea- 
surement method for electronic states [Q. We set e—1 
and fee = 1. 

The fundamental idea of the quantum-dot qubits based 
on charged states is as follows. A qubit is composed of 
two quantum dots coupled via a thin tunneling barrier 
and a gate electrode that is attached on a thick insulating 
material. This is a bistable well structure [|9| where the 
electronic quantum state of the coupled dots is controlled 
by the gate bias. First, we consider the one-dimensional 
arrayed qubits(Fig. y). The quantum dots are, e.g., Si 
nanocrystals M-E3] or GaAs dots |0,n. The quantum 
dots are assumed to be sufficiently small for charging 
effects to be observed. N^i and Npi are the numbers 
of excess electrons from the neutral states in two quan- 
tum dots. One excess charge is assumed to be inserted 
from a substrate first and to stay in the two-coupled dots 
{Csi > Cci). When the excess charge exists in the up- 
per dot and lower dot, we call them |0) = 1 1) state and 
|1) — II) state, respectively. The quantum logic gates are 
assumed to be operated in the region where the |0) and 
|1) states are near-degenerate, then the system becomes 
a two-state system raj|,p3| . The qubits are arranged one- 
dimensionally and are capacitively coupled. 

Here we show 

the Hamiltonian of the one-dimensionally arrayed cou- 
pled quantum dots for quantum computation by starting 
from the tunneling Hamiltonian: 
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where Oj (bi) describes the annihilation operator when 
the excess electron exists in the upper (lower) dot, and 
fai (£/3i) shows the electronic energy of the upper (lower) 
dot. There is no restriction on the number of coupled 
quantum dots, N. Hch is the charging energy that in- 
cludes the interaction between qubits. Because we con- 
sider the Coulomb blockade in the weak coupling region, 



the resistance of the interdot tunnehng barrier should be 
larger than Rk{= h/e^ ^ 25.8kJ7) and the operational 
temperature be less than the charging energies. Thus, 
the operational speed should be less than the CR con- 
stant of the capacitance network so that the double well 
potential profile generated by the charging energy is effec- 
tive (adiabatic region) . The criterion of the operation is 
discussed below in detail. As for the interaction between 
qubits, the distribution of the excess charge is considered 
to be antiferromagnetic due to the repulsive Coulomb in- 
teraction. We can show that this interaction between 
qubits is an Ising interaction by minimizing the general 
formulation of the charging energy: 



where \]i)^\n.i = 1) and | |j) = \ni = -1), U^ = {a\a^ 
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with constraints -Nai = qM-qm+qiyi-qiyi-i+qE.i-qFi-i+ 
qm-i+qii and -Niii = q^i-qci+qGi-qGi-i-qEi-i+qFi 
{i — 1,..,N). The last line of Eq.(||) shows the effects 
of other gate electrodes on the ith qubit (cross-talk) . By 
using Lagrange multiplier constants, the energy of the 
system is obtained as a function of the relative excess 
charge rii = Na^ — Nj^. and the total excess charge of the 
two quantum dots Ni = Na^ + Ni3.{=l ): 
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with Cat = Caz + Ccz + 4Cbz +2{CDi + CEi) + Cm + 

Cli, Cbi = CAi+Cci+'2{CDi+CEi) +2{CDi-l-\-CEi-l) + 

Chi+Ch, Cci = Cci— CAi+CHi+ Cli, Cdi = Cm+CEi, 
Cei = Coi — CEi {Coi = Cci, Csi = Cpi) and Di = 
CaiCbi — C^i and Qv = CAiVgi + Cni-iVgi-i+CjiVgi+i. 
It is assumed that the coupling between qubits is smaller 
than that within a qubit and we neglect higher order 
terms than (Cjj/Di)^(<C 1). This assumption is valid 
when the distance between the quantum dots in different 
qubits is larger than that between quantum dots in a 
qubit. From Eq.(|^), we define the characteristic charging 
energy of the system as Ec = Cb/{2D). We consider the 
gate voltage region, where the ni=-l state and ni = l state 
are near-degenerate as in Ref. [0j8||2^ , and we obtain the 
Hamiltonian of the coupled quantum dot system: 
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- [CbiCbi+lCei+1 + CciCci+lCdi+l] ■ (6) 



Qy^ includes e^i and e^i and shows the gate voltage when 
the |0) and |1) degenerate (on resonance). We control the 
time-dependent quantum states of the qubits in the vicin- 
ity of on resonant gate bias by applying a gate voltage 
such as Vi (r) = V^''^^ + Vi (r) where V^°^ is the gate bias 
of on resonance. In this on resonant region, a transfor- 
mation of the coordinate: 
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is convenient (we neglect the phase of U for simplicity). 
Then the Hamiltonian Eq.(0) can be described as 
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where 
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Svicrir) = [Cm-iVi-i{T)+CiiVt+i{T)]/CAi is a cross-talk 
term and is an effect of other gate electrodes. Hamilto- 
nian (H) is an NMR Hamiltonian when we regard uji = 2ti 
as a Zeeman energy and Ai(r) as a transverse magnetic 
field, if w > Ao ^ J. Thus, when Ai(r) = Aoi cos(wir -t- 
5i) and we take a rotating wave approximation by the 

unitary transformation C/rwa(''') ~ cxp(— i ^.-^-^ cj^r/j'^), 
we have 

^ Aio -^ 
Hrwa= - E^ K» COS 5, + /^ sin 5,] 
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The pulse process is carried out when the oscillating elec- 
tric field is applied (Aq / 0) where the interaction term 
is considered to be able to be neglected because Aq ^ J. 
For example, Tr^-pulse in the (Si, bi) basis corresponds to 
TTaj-pulse representation in this (a+i, a_i) basis and is car- 
ried out when i5 = and Aor/2 = tt. In general, rotation 
Rj{6) = e '^^"' (7 = x,y,z) in the {ai,bi) basis is inter- 
preted by R[{9) — UQ^Ri{9)Uo in the {a+i,a-i) basis. 



The evolution process is carried out when the oscillating 
voltage is not appHed (Aq =0). A two-qubit CNOT gate 
is described by a pulse sequence of the form t/^j^oT °^ 

t/o-^i?^, (7r/2)i?;.,y (7r/2)i?;, (V2)i?., (-7r)i?;, (-7r/2) t/o 



The 2-body interaction _Ry(— tt) 



-^/../.4 



is obtained by using an average Hamiltonian theory p7| . 
Here we apply the "Carr-Purcell" sequence t — tTj, — 2t — 
TTy — T to the Hamiltonian Eq. (@) in order to average the 
effect of Zeeman term to zero. Thus, we can show that 
the Hamiltonian of the coupled quantum dot system can 
be reduced to that of the weakly coupled system in NMR 
and quantum operations can be carried out in a manner 
similar to that of operations in NMR quantum comput- 
ers. 

Similarly, we can show that the Hamiltonian of the 
two-dimensionally arrayed coupled dots is reduced to 
that of an artificial Ising system controllable by elec- 
tric fields. Here we consider the case where there are 
only nearest-neighbor capacitive couplings {Ce and Cp 
are neglected in the above formulation). If we express 
the coupling strength between the sites i = (i,j) and 
i+x = {i + l,j) as J^ and that between the sites i and 
i+y = (^,J■ + 1) as Jf (1 < « < N,, I < j < N_,: 
N = N^Ny), we have a Hamiltonian similar to Eq. (H) 
where 
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- (CbiCbi+y + C'ciCci+y)C' jj-^, (13) 
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with A EE CaiCbi - C^i, Cai = C^i + Cci + 4CBi, Cbi 

= CAi + Cci, Cai = Cci-CM- CSi(= Cgi) and Cl,{= 
Cgj) are capacitances between qubits in x-direction and 
y-direction, respectively. 

Here we summarize the criterion for realizing the 
above-mentioned scheme. First of all, the Coulomb 
blockade should be effective at the operational tempera- 
ture T and we have T <^ Eq- In view of time-dependent 
operation, the excess charge is assumed to be affected 
by the electronic potential generated by the capacitance 
network. Therefore, all quantities concerning time evolu- 
tion should be smaller than the CR constant of the net- 
work. Here we take CR = Cint^int where Cint(= DjCb) 
and _Ri,it are capacitance and resistance of the interdot 
tunneling barrier in a qubit, respectively. By including 
the condition for using the pulse sequence as mentioned 
above, we have a condition for the operation: 



T < J < Ao < t < {CR)- 



(14) 



We can roughly estimate this criterion by taking typi- 
cal values, vq =2.5 nm (radius of a quantum dot), dyi=8 
nm, dB= 1.5 nm, dc= 2.5 nm, and the distance between 
qubits du is 12 nm (eox =4 (Si02) and esi=12), reflect- 
ing several experimental data ||20|-|2^. Using relations, 
CA,c=27reox'ro/(rfA,c+(eox/esi)?'o) and Cs = 27reox?'o/(dB 



0.4meV and J~0.1 meV. If i?int is of the order of M17, 
we obtain (Cint-Rint)"^ ~ S.lTHz and greater than t ~ 
lOOGHz. Thus, the condition ( p^ ) is satisfied if the oper- 
ational temperature should be much less than IK. If we 
could prepare r^ — 0.5nm, dB = 1.2nm and do = 2nm, 
we obtain t ^^120 K and J ~ 90K and quantum cal- 
culations can be expected to be carried out at around 
liquid nitrogen temperature. The effects of cross-talk are 
of the order of Cui/CAi ~ dAi/doi and cannot be ne- 
glected even when gate electrodes are set closer to the 
corresponding quantum dots. However, we can control 
the cross-talk effects by adjusting Vi{T) to obtain the re- 
quired Ai(r). 

At the near-degeneracy point, the system becomes a 
two-state system and the estimation of the decoherence 
(in order of /isec) discussed in Ref. 1l2] may be applicable. 
Until now, we have not known of any corresponding ex- 
perimental data for the decoherence time [M . The point 
is that various methods developed in NMR spectroscopy, 
such as the composite pulse method [0, can be utilized 
to reduce the imperfections of the pulse and coherence 
transfer, that is, errors that are brought about in the 
operations. In addition, if the speed of quantum compu- 
tations can be increased to much more than the shortest 



decoherence time, t^. 



^{^ 10 ^"^ s), group-theoretic 



approaches [^^ for decreasing the decoherence will be 
effective. In this case, we have to reduce the correspond- 
ing CR of the junction such that Eq.([l4|) holds and the 
cycle time of operations should be much less than Tc- 
Thus, the errors and the effects of decoherence of a cou- 
pled quantum-dot system will be reduced by developing 
many contrivances. 

Hereafter we consider the one-dimensionally arrayed 
qubits for simplicity. Next, we quantitatively illustrate 
the reading out process based on an FET structure 
(Fig.|l|). Measurement is carried out, after quantum cal- 
culations, by applying a finite bias Vd between the source 
and drain |12|. The detection mechanism is such that 
the change of the charge distribution in a qubit induces 
a threshold voltage shift AVth of the gate voltage above 
which the channel current flows in the substrate. AVth 
is of the order of edq/coy^idq is a distance between the 
centers of two quantum dots in a qubit) ||21|]. The effect 
of AVth differs depending on the position of the qubit, 
because the width of the depletion region in the sub- 
strate changes gradually from source to drain. A sim- 
plified model of a metal-oxide semiconductor field-effect- 
transistor (MOSFET) including velocity saturation ef- 
fects (9 in Eq. (p^)) pq is used for the current that 
flows under the ith qubit 
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-f-2(eox/esi)'^o), we obtain Ec 



13meV(150K), t 



where A = ZiiqCq/Lq (Z is the channel width, /iq is 
the mobility, Lq is the channel length of one qubit, and 
Co is the capacitance of the gate insulator) , r]i = 1 + Q 
where Ci is determined by the charge of the surface de- 



pletion region, and Vi is the voltage of ith qubit to be 



determined from V^ 



Vb and /^^^ 
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The detection should be carried out before the change of 
the quantum states of the qubits. For this purpose, A 
should be as large as possible |Q . The threshold voltage 
of ith qubit Vthi is given by Vthi = Vth(Cj) + ^^thj- Fig- 
ure 0(a) shows the ratio of the current change |/i-/o|//o 
as a function of Vd in 8 qubits with Q = 0, where Iq 
is the initial current and li is the current when Vgi-Vthi 
of ith qubit changes by 10%. This ratio is largest for a 
qubit near the drain (z = 8) because it has the narrow- 
est inversion layer. To show how to distinguish qubits, 
we compare the current where only the ith qubit shifts 
its threshold voltage with that where only the (i-l-l)th 
qubit shifts its threshold voltage. Figure ||(b) shows the 
results for (i)i=l, (ii)i=A^/2, and (iii)i=7V — 1, with the 
same voltage shift as in Fig.g(a). The maximum allow- 
able number of arrayed qubits depends on the sensitivity 
of the external circuit to the channel current. When the 
density of acceptor in the substrate is of the order of 
10"'^^ cm~'^, the number of acceptors below one qubit is 
less than one. If this effect is represented by Q as a ran- 
dom number, the above ratios may increase or decrease 
and the overall features are similar to those in Fig.|| To 
construct a large qubit array, additional dummy qubits 
are needed over the source and drain so that separated 
qubits on different FETs are connected. 

In conclusion, we have theoretically shown that the 
Hamiltonian for a weakly coupled quantum-dot array in 
the Coulomb blockade regime is reduced to that for NMR 
spectroscopy. The flexible quantum information process- 
ing developed in the NMR quantum computer and a va- 
riety of error-correction methods in conventional NMR 
are applicable to the quantum-dot system that is con- 
sidered to be the most feasible system in view of the 
present technology. The difhculty of scalability in the 
NMR computer, namely that the overlap of pulses re- 
stricts the number of qubits H], is overcome in the op- 
erations by individual gate electrodes. The disadvantage 
of the short decohcrcnce time is also compensated for by 
another advantage, namely that the measurement proce- 
dure is compatible with classical circuits. It is expected 
that a transmitting loss of signals to classical circuits will 
decrease due to this compatibility. Recently, nanofab- 
rication of two-dimensionally distributed self-aligned Si 
doubly- stacked dots has been successfully realized in the 
form of a non-volatile memory device p2[ and the de- 
tailed analysis of the behavior of electrons, such as an 
artificial antiferromagnet, is expected to be performed. 
The iV-qubit system composed of arrayed quantum dots 
is expected to be developed as a result of further advances 
in fabrication technologies. 

The author thanks N. Gemma, R. Katoh, K. Ichimura, 
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FIG. 1. Schematic of the capacitance network of a 
one-dimensionally coupled quantum-dot array. The capaci- 
tances Caz, Cbi, and Cci, the gate Vgi and the two quantum 
dots constitute the ith qubit. The excess electron moves be- 
tween the two quantum dots in a qubit via the tunneling 
barrier, and electron transfer between qubits is prohibited. 



FIG. 2. /; represents the current where only the ith 
qubit shifts its threshold voltage Kh(0)+AKhi- lo represents 
the current where all qubits have the same threshold voltage 
Vth(O). (a) The ratio of current change [/; — Io\/Io as a func- 
tion of Vd in an 8-qubit quantum computer (i = l, i = 4, and 
i — 8). (b) The ratio of change \Ii — Ii+i\/Io as a function of 
the number of qubits, in the case i = l(near source), i = N/2 
(middle) and i = N — l(near drain) at Vd ~ 1.5V. The ef- 
ficiency of detecting the quantum state of qubits is highest 
when the qubits near the drain change their states. In (a) 
and (b), Vg - Kh=2V and 9=0.3 V"^ in Eq.(7) in text. The 
threshold shift is 10% of Vg - Kh- 
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